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Abstract 

We examine the Onsager algebra symmetry of r^-matrices in the superintegrable chiral 
Potts model. The comparison of Onsager algebra symmetry of the chiral Potts model with 
the s/2-loop algebra symmetry of six-vertex model at roots of unity is made from the aspect of 
functional relations using the Q-operator and fusion matrices. The discussion of Bethe ansatz 
for both models is conducted in a uniform manner through the evaluation parameters of their 
symmetry algebras. 

2003 PACS: 05.50.+q, 02.20.Uw, 75.10Jm 
2000 MSC: 14H70, 17B67, 82B23 

Key words: Chiral Potts model, Six-vertex model, Onsager Algebra, sfe-loop algebra 



1 Introduction 

The symmetry of quantum spin chains and the related lattice models has recently attracted cer- 
tain attention due to their close connection with diverse areas of physics as well as mathematics. 
However, up to the present stage, only limited knowledge is available about the symmetry of lattice 
vertex models, and few exact results are obtained in this area. Even the s/2-loop algebra symmetry 
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of the six- vertex model at roots of unity, found in [11], has not been fully understood till now, given 
that much accomplishment has been made on the study of evaluation parameters related to the 
symmetry algebra representation. Some conjectures supported by numerical evidence remain to 
be answered, (see [12] [14] and references therein). Though the understanding of the symmetry of 
eight-vertex model in [15] [16] is still rudimentary in the present stage, the conjectural functional- 
relation-analogy discovered in the study on the eight-vertex model and the TV-state chiral Potts 
model (CMP) 3 in [7] did lead to exact results about the Onsager algebra symmetry of r^-models 
in the superintegrable CPM [20] [21]. In the study of CPM, the lack difference-property for the 
rapidities is considered as the characteristic nature which distinguishes CPM from other known 
solvable lattice models. Nevertheless, progress made on the transfer matrix of CPM for the past 
fifteen years, culminating in the recent Baxter's proof of the order parameter [6], has provided the 
sufficient knowledge for the understanding of the detailed nature about the symmetry of superinte- 
grable CPM. By this, the study of CPM could suggest a promising method to help the symmetry 
study about the six-vertex model at roots of unity as the limit case of eight- vertex model from the 
approach of functional relations, a scheme proposed in [16]. In this paper, we examine the similar- 
ity of the symmetry structure of two lattice models: the superintegrable CPM and the six-vertex 
model at roots of unity. The symmetry of superintegrable CPM is described by the Onsager alge- 
bra, obtained in [20] [21], with a short explanation in Sec. 2. The six- vertex model at roots of unity 
possesses the s/2-loop algebra symmetry by the works in [11] [12] [14]. We present a discussion of 
symmetry of six-vertex model, parallel to the theory of CPM, through Bethe roots and evaluation 
parameters in Sec. 3, and give some concluding remarks in Sec. 4. 

Notations. In this paper, Z, C denote the ring of integers, complex numbers respectively, N 
is a positive integer > 2, Z^v = Z/iVZ, and i = y — 1, oj = e n . Denote by X, Z the Weyl operators 

of the vector space C^, defined by X\n) = |n + 1), Z\n) = uj n \n) for n € Zjv, and <8> C N is the 
n-tensor product of C^. 



2 The TV-state Chiral Potts Model 

2.1 Rapidity and functional Relation of chiral Potts model 

In the study of CPM as a descendent of the six-vertex model, Bazhanov and Stroganov obtained 
the following 3-parameter family of Yang-Baxter solutions for the inhomogeneous R-matrix of six- 
vertex model, 

/ tu> - 1 \ 

t- 1 uj - 1 

t(u-l) uj{t-l) 
\ tu - 1 J 



R(t) = 



with the C -operator entries parametrized by a four- vector ratio p = [a, b, c, d] [8, 20]: 



b 2 G p (t) = 



b 2 - td 2 X {be - ujadX)Z 
-t(bc - adX)Z- 1 -tc 2 + ua 2 X 



t € C , 



which satisfy the following Yang-Baxter equation: 

R{t/t'){G p {t) (g) 1)(1 (9) G p {t')) = (10 G p (t'))(G p (t) l)R(t/t f ). 



(1) 



3 In this paper, the discussion of the chiral Potts model will be confined only to the full homogeneous lattice by 
taking p — p' in [7] . 
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Hence the same relation holds for the monodromy matrix of size L, ®f =1 G Pt e(t) with G Pt e(t) = 
G p (t) at the site I. Therefore, the Tp^-matrix, 

L 

ij 2) (^=tr««®^(wt)) for tGC, (2) 
l=i 

form a commuting family of (<g C^-operators. The Zjy-operators X, Z of C N at the site j give rise 

to the Weyl operators XjZj of (g with the relations: ZjXj = uA^XjZj, [Zj, Z,-] = [Xj, Xj] = 

and Z^ = Xj^ = 1. Then the spin-shift operator of (g C^, X := Ylf =1 Xj, defines the Z^-charge 

Q, and commutes with r p 2 \t). 

The rapidities of iV-state CMP are elements in the genus (iV 3 — 2N 2 + 1) curve 2U in the 
projective 3-space P 3 , defined by the equivalent sets of equations: 

„ / ka N + k'c N = d N , _ f a N + k'b N = kd N , f kx N = 1 - Vfi- N , 

W: \kb N + k'd N = c N ^ \ k'a N + b N = kc N < ^\ky N = l- k'n N [6) 

where [a,b,c,d] £ P 3 , (x,y,fi) = G C 3 , k' is a parameter with k 2 = 1 - k' 2 / 0,1. By 

eliminating ^ in the last set of above equations , and using the variables t := xy,X := fi N , one 
arrives the hyperelliptic curve of genus N — 1, t N = — - ~^^T fc A — as a X 2 -unramified quotient 
of (3). The rapidities possess a large symmetry group, in which the following two will be needed 
in our later discussion, 

U : [a, b, c, d] i-> [wo, 6, c, d], (x, y, /i) i-> (wx, y, //), (t, A) i-> (wt, A); 

C : [a,b,c,d] ^ [b,a,d,c], (x, y, fi) i-> (y, x, /i" 1 ), (t, A) i-> (t, A -1 ). 

The Boltzmann weights Wp i? , of the CPM, depending on two rapidities p,q £ W, are 
defined by 

Wp, 9 (ra) _ T-rn d v b q -a p c q u)i W p , q (n) _ yin ua v d q —d v a q ui 
Wp'qiO) ~ Hi=l b p d q -c p a q u)i ' W p ,,(0) ~~ 11 -5 ,=1 Cpb q -b p c q u3 ' 

which satisfy the A^-periodicity property for n by the constraint (3) of rapidities. The CPM transfer 
matrix of size L with periodic boundary condition, L + 1 = 1, is the ((g) C Ar )-operator defined by 

L 

T cp (p;q)aC:,* L L = \{W M (a l -a' l )W p ^a l -a' l+l ) , a l ,a' l eZ N . (5) 
z=i 

For a fixed p € 2U, {T cp (p;q)} qG %sj form a commuting family of operators by the following well- 
known star-triangle relation of Boltzmann weights: 

N-l 

]T W qr (a' - n)W pr (a - n)W pq {n - a") = R pqr W pq (a - a')W pr {a' - a")W qr (a - a"), 

n=0 



where R vqr = M^l w ith f vg = ( ^M^WiJ)) ) . Then T c Jp;q) commutes with X and the 

/pr v n„= w Pi( n ) j 

spatial translation operator Sr (which defines the total momentum P £ Z^). Denote T cp (p;q) := 

(2) 

T cp (p; q)SR. The transfer matrix T cp (p; q) can be derived from t p (t q ) with p G 22} as the auxiliary 
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"Q" -operator, as discussed in [2] on the TQ-relation of the eight-vertex model. One arrives T^T cp - 
relation ((4.20) in [7]) 4 using the automorphism U in (4): 

ri 2 \t q )T cp (p;Uq) = Vp {q)T cp {p{-q)+i ?p {Uq)XT cp {p ] U 2 q) (6) 
4 2 \t q )= ^(?)T cp (p; g ) + ^([/<|)ir cp (p; [/ 2 g)j T cp (p; [Z?)- 1 (7) 

where <p p (q) := ( ^ggKffi) )^, ^ p ( g ) : = (^|g^£zM)^. By (7) and the commutativity of 

T cp (p;*), [rp 2 ^ (i g ), T cp (p; </)] = for p,q,q' £ 2U. The fusion operators T^\t) for < j < N are 
determined by the following T cp T cp -relation for < j < N, ((3.46) for (l,k) = (j,0) in [7]) with 

rf := 0, r« := /, 

T cp (p; q)f cp (p; CWq) = r p , qhj ., p , q (f (f ff ) + ^ -j) ^ (g) 

where z{t) - y ^ J , a p {\ q ) - y Xq(1 _ k , Xp) , ) and r p , q - [ (x n _ x N )(y N _ y N )(tp _ tq) ) , 

h r,P,g = (it 1 ! (y p -L™l~)"t v -i™t q ) ) ■ B y ( 8 )> one can derive the fusion relations of r^s ((4.27) of 
[7] ): 

4 J \t)4 2 \u^H) = z^-^Xrt^it) + 4 3+1 \t), l<j<N, 

4 N ^\t) := z^Xr^^t) +u(t)I { ) 

where u(t) := a p (X) + a p (A _1 ). Note that with r p 2 \t) in (2) for p £ P 3 , the validity of fusion 
relation (9) provides a characterization of the rapidity constraint (3) for p £ 2U, (Theorem 1 in 
[20]). Using (7) and (9), one can express in terms of T cp (p;q), hence the r^^T cp -relations for 
1 <j <N + 1, ((4.34) in [7]): 

4 j \q) = T cp (p; q)T cp (p; Wq) E J m =o U P {q)<Pp{Uq) ■ ■ ■ M^'M* 

V x (10) 

<^(£/ m+1 <z) • • ■^(^- 1 g)T cp (p; U™q)- l T cp (jp- U^q^X^^J . 

By substituting (10) in (8), one obtains a single functional equation of T cp (p;g) ((4.40) of [7]): 

?c P Cp; g) = E™=o C m - p (q)T cp (p; q)T cp (p; U m q)- 1 T cp (p- U m+1 q)' 1 X~ m ~\ (11) 

where C m , p {q) = <p p (q) ■ ■ ■ ^(U^q^U^q) ■ • • <p p {U N ^q) ( N %^^%) q) ) L '• 

2.2 Bethe Ansatz and Onsager algebra symmetry in superintegrable chiral 
Potts model 

For CPM in the superintegrable case, i.e., the rapidity p given by [x p = 1, x p = y p = 772, where 
rj := (j^r) 7r ) simplification occurs for the functional relations. We shall omit the label p appeared 
in all operators for the superintegrable case, simply write T^(t), T cp (q), etc. As q tends to p, to 
the first order of small e, one has ((1.11) in [1]), 

T cp (q) = 1[1 + e(JV -1)L] +£#(*') 
4 Tp 2 \q),T cp (p;q) in this paper are the operators r^ , T g in [7] respectively. 
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where H(k') = Hq + k'H\ is the Zjv-quantum chain in [17, 18], with the expression: Hq = 

-2EiiEn=i 1 T=2=n, Hi = -2EiiEl"i 1 which satisf y the Dolan-Grady relation [13], 

hence generate the Onsager algebra representation where only the spin-i subrepresentation oc- 
curs as irreducible factors [1, 3, 5, 10, 19]. Through the gauge transform by M = dia. [1,772], the 
monodromy matrix G(t) in (1) becomes 

G(t)-( ~ l ~~ tX ^Z^ 2 ) t-n-H 

which is again a Yang-Baxter solution (1) for a = b = c = d = 1. Hence T^ 2 \t) = T^ 2 \t), where 
r( 2 )(t) is the trace of the L-size monodromy matrix associated to G(ut). Write T^\t) = f^\t), 
the fusion relation (9) has the form: 



T^\t)f^((J-H) = (1 - (J~H) 2L tV-V® u l X + r^ +l \t), 1 < j < N , 

f (N+l) ^ = (1 _ £)2L ~(JV-1) (^) W L X + 2(1 — t N ) L . 



(12) 



By examining commutators of Hk with the entries of the monodromy matrix constructed from 
G(t), one can show [Hk, (t)] = for k = 0, 1. It follows the Onsager algebra symmetry of re- 
model (Theorem 1 of [21]). However, the understanding of the detailed nature of Onsager algebra 
symmetry in the superintegrable CPM still requires the full knowledge about eigenvalues of CPM 
transfer matrix, which was solved by the Bethe-ansatz method in [1, 4, 5] as follows. 

For parameters v\, . . . ,v mp with (—Vi) N / 0, 1 and Vivj 1 7^ 1, u for i 7^ j, consider the rational 
function 

V(t) = uo- p > Y ( I ( 1 F(t) := TT(1 + uvit) (13) 

where P a , Pb are integers satisfying < r(:= P a + Pb) < N — 1, Pb — P a = Q + L (mod N). V(t) 
is invariant under t 1— > wi, hence depending only on t N . The criterion of V(t) as a t-polynomial is 
the following constraint for VjS, ((4.4) in [1], (6.22) in [4]): 

- _ 2 ) = 11 . » = l,...,Tnp. 14 

TJj + W ^ - UWj 

Here the non-negative integer m p satisfies the relation LPb = m p (Q — 2Pb — m p ) (mod N). The 
total momentum P is given by e lP = uj~ Pb IlSi I+^J~- The above relation is indeed the Bethe 
equation of r^-model (Theorem 3 of [21]). Then V(t) is a simple t^-polynomial of degree niE = 
y (N-i)L-r-2m p ^ w ^ negative real roots, (Theorem 2 of [21]). Let s\,...s mE be the t^-zeros of 

V{t), and define G{\) = Upi x+1±< £r 1)w i where Wj := ( Sj ) ^ • Then = ^(A)g(A- 1 ). One 
has the following expression of T cp (g)-eigenvalues ((1.11) in [1] and (21) in [5]): 



T CP (1) = N L ^ N Xq N ^ {^x q )^{^y q f^- p ^G{X q ) (= e^fo)), (15) 

(rj— -\) L ^i 1 ) 

which gives the energy value of H(k') ((2.23) of [1]): 

E = 2P^ + Nm E - (TV - 1)L + k' UN - 1)L - 2P M - iVm £ + 2(P b - P a ) \ + N(l - k') ^ ± Wj . 

V / j=i 
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Therefore the t( 2 ) -degenerate states associated to the Bethe roots V{S form an irreducible Onsager- 
algebra-representation space of dimension 2 niE , which we associate the following normalized CPM 
transfer matrix 5 : 

m _ w (=: e -.pQ (5)) . 

N L (1 - Tj 2 x q ) L (r] 2 x q ) Pa (r] 2 y q ) p b[i q M 
The Q-eigenvalues and the functional equation (11) now become 

0(9) = T§ G ^ ^ C<? ) = eiP T§ G ^ ( 16 ) 

Q(Cq) = # Em=J ^g^Q(c?)g(^ m g )- 1 Q(f/ m + 1 g)- 1 . (17) 

The relations, (6) (10) (8), now become 

T<V(t q )Q(Uq) =uo- p «(l-t q ) L Q(q)+u p »(l-Lot q ) L Q(U 2 q), (18) 
rWfo) = W 0--D^ E ^ ^ II^Z^ Q (g) Q ( ^ g) -i Q( ^ g)Q( ^ + i g) -i^-^ j (19) 

= ^gg^ rfa-)^,) + ^g^L r(^)(^T g ), < j < N. (20) 

1 lm=0 ^ 1 lm— j ^ U 

By (16), the relation (18) yield ((6.18) in [4]): 

T^(t q )F{Jt q ) = u- Pa {\ -t q ) L F{t q ) + uj Pb {\ - uji q ) L F{oj 2 t q ). (21) 
Using (19) and (20), one obtains the f^-polynomial expression (Theorem 3 (ii) of [21]): 

t^)-^ u H fc=0 (l - w L m=0 ' 2<j<iV, 

?F(i)F(^7)P(i) = ""-^-^ rO-)^ + "T-l^-M ^^t), < j < N. ^ 

3 Six-vertex Model at Roots of Unity 

3.1 Bethe equation of the six- vertex model 

L „ 

The transfer matrix of the six-vertex model of an even size L is the (<8> C )-operator constructed 
from the Yang-Baxter solution 6 

-<T Z 

q 2 — 2 s 

(q — q _1 )<T + , z^q~ — z~q^ 



Z 2q 2 - z 2 q a 2 (q - q x )cr_ . , 

L = z _, _ z , o ,cr± : Pauh matrix, 



for the R-matrix 

/ z^-q-z^q- 1 

-1 1 

z 2 — z2 q — q 

— 1 =1 1 
q — q z 2 — z2 



J2(z) 



z^"q - z^q- 1 ) 



5 The Q-operator here differs from the Q cp in [21] by a scale factor: x q ) Pa (r/~ y q ) Pb fj, q Ptl Q(q) = Q cp (q)- 
6 Here we use the convention in [14], Eq. (1.3): a = isinh — i7), b = — isinh |(v + i7), c = — isinhi7, with the 

variables z = — e v ,q = — e 17 . Note that the q here differs from the q in [14] by minus sign for its connection with 

U q (sh) as Eq. (2.3) in [12], where — e~ z is equal to 22 here. 
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as the trace of monodromy matrix: T(z) = tr aua; (0^ =1 Lg(z)) for z G C. The logarithmic z-^- 
derivative of T{z) at z = q _1 gives rise to the XXZ chain with the periodic boundary condition: 

1 L 1 

Hxxz = - » £fo 1<7 i+i + ^f+i + A<r|af +1 ), A = -(q + q x ), 

which has been a well-studied Hamiltonian initiated by Bethe in 1931 [9]. A major result of the 
study is that the ground state energy for the value S z {= | J2e a f) ls determined by an appropriate 
solution of the following Bethe equation for the variable v = —z~ l : 

r» +q ~V = -q^f r~ q ~ 2 ^ rn=±-\S Z \. (23) 

The Bethe-equation technique was further extended to the method of Baxter's TQ-relation in eight- 
vertex model; when applying to the six-vertex model, there exists a non-degenerated commuting 
family of Q-operators with the following relation (see, Chapter 9 of [2]): 

T(z)Q(s) = q~ 2 l sZ l(l - qz^QiU'h) + (1 - q- l z) L Q(Us). (24) 

Here s is a suitable multi- valued complex coordinate of z, and U is a s-automorphism which induces 
the z-transform: z i— > q 2 z. Note that there are many such Q-operators, however all give the same 
Bethe equation (23) through Eq.(24). 

3.2 Evaluation polynomial and fusion relation of six-vertex model at roots of 
unity 

For the root of unity case with q 2JV = 1, i.e. q 2 = u>, the six-vertex model possesses the s/2-loop 
algebra symmetry [11]. The Bethe state corresponding to the Bethe roots is the "highest weight" 
vector of an irreducible representation of the s/2-loop algebra, with the evaluation parameters char- 
acterized by the Drinfeld polynomial [12, 14]. By studying the creation s^Toop current operator in 
the ABCD-algebra, the Drinfeld polynomial for a Bethe root {vj}^ of (23) is given by Eq. (3.9) 
in [14]. Denote t = qz, and define the integer r by r = ^ — m (mod N), < r < N — 1. The 
Drinfeld polynomial is indeed the t^-polynomial associated the following polynomial P(t) 7 , 

P(t) = v LgMr , F(t) :=f[(l + q~\t), (25) 
which has a similar form as (13). Indeed with F(t) in (13) or (25), let H(t) = "y^p 1 — t in 



CPM, six-vertex model respectively. The function P(t) defined by P(t) = J^jLo 1 ^7^p(^7+i7) ^ s 
invariant under t i-> ut. The condition on roots of F(t) so that P(t) is a polynomial is provided 
by Bethe equation (14), (23) respectively. Define the T( 2 )-operator by T^ 2 \t) = z^T(z) in the 
six- vertex model, and - bi - 1 ^ t ' 1 t( ^Jf ^ in CPM. Then Eqs. (18), (24) are combined into one 
T^Q-relation: 

T^(t)Q(q) = W - r F(t) L Q([/- 1 g) + H^-^QiUq), U N = 1. (26) 

In [14], the evaluation function is given by (3.9) there: Y(v) = E™=o + 2 U + Vh), wher e a(v) = 

in (2.47) of [14]. In terms of variables z,t and Bethe roots vis here, a(v) = 



J~J sinh ■i- (vi — v) sinh ^ (v^ — v-\-2i-y') 

^- L (UZi v ^~ lz )^ +m YY^S&^y which implies Y{v) = 22m_L (nr =1 « i )*^ +m+r i'(*) 
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The -operators are defined recursively through the following fusion relation for j > 1 by setting 
T (0) = Qj T (i) = H (u- i: t) L , 

r 0)(t) T ( 2 )(^- 1 f) = uj-'HiJ-^T^-^it) + H{uj^ 2 t) L T^ +1 \t). (27) 
By (26), the induction-argument yields the T^Q-relation for j > 0: 

j-i 

T&(t) = Q(U' 1 q)Q(U^ 1 q)J2(^ kr H(io k - 1 t) L Q(U k - 1 q)- 1 Q(U k q)- 1 ). (28) 

fc=0 

By this, one obtains the boundary condition of the fusion relation: 

r (*+i)(t) = u^T^^iujt) + 2H(u~ 1 t) L . (29) 

In CPM case, with the identification T^(t) = ^'^(i-^)^^' 1 *) E (27)-(29) are the 

same as Eqs. (12), (19). While in six- vertex model, the fusion relation and T^Q-relation hold for 
any Q-operator satisfying T( 2 )Q-relation (26). 

For a polynomial F(t) with Bethe roots v jS, by (28) the corresponding T^> -eigenvalue is deter- 
mined by the relation 

T^{t)F{t) = u^Hit^Fiuj-H) + Hiuj-^Fiui). (30) 

Using Eq.(27) to express T^ +l ^ in and T^ -1 ), by induction argument, one obtains the form 
of rC?) -eigenvalues from Eq.(30): 

3— 1 rr/ ,fe— 1T\L — fcr 

^^F( W "WiEM^> J>1, (31) 

which implies 

t'Fiuj-^FiJ'^Pit) = T ij \t) + w^'T^-^O^t), < j < N. (32) 

Eqs.(30)-(32) in CPM case are the same as Eqs. (21), (22). Note that Eq.(32) is a consequence of 
the QQ-relation (20) in CPM case, which encodes the detailed nature of Onsager algebra symmetry 
in the derivation of Eq.(16). However in the case of roots-of- unity six- vertex model, the QQ-relation 
has yet been found, even though the s^-loop algebra symmetry together with evaluation parameters 
has already been known [11, 12, 14]. Based on the understanding in the CPM case, we now describe 
a similar, but speculated, structure about the QQ-relation in six-vertex case as follows. Consider 
the curve W and its symmetries U, C: 

W:w 2 = t N , U : (w,t) i-> (w,ut), C : (w,t) i-> (-w,t). 



For odd N, the above curve is parametrized by s = ia, and the automorphism tp : s i— ► qs 
gives rise to the above symmetries by U = tp~ 2 ^~\ C = (p N . The polynomial P(t) in (25) are 
expected 8 to have the simple f^-roots s\, ...,sm 

with P(0) + 0. Define G[w) = U.jLi(y/sJ ~ w), 

then G ^ V Q^i W ^ = 7J§)- ^ n tne eigen-space of T(z) corresponding to F(t) determined by a Bethe 
root, The Q-operator has the following Q-eigenvalues: 

Q(9) = F(f) ^' Q{Cq) = F{ ^^GW iorq = {wJ) - 

The above conditions reveal the s^-loop symmetry of six- vertex model, as the role of Eq.(16) for 
the Onsager algebra symmetry in CPM. Hence such a Q-matrix, if exits, must possess certain 
constraints in order to incorporate the symmetry of six-vertex model as discussed in [2] Sees. 
9.1-9.5. 



8 The statement is true in the case r = 0, which has been justified in [12]. 
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4 Discussion 

In this paper we have examined the symmetry structure of the superintegrable CPM and the six- 
vertex model at roots of unity by the method of functional relations. For the superintegrable 
CPM, exact results about the Onsager algebra symmetry of the r^-models are obtained using 
the explicit form of eigenvalues of the CPM-transfer matrix. Based on common features related to 
evaluation parameters of the symmetry algebra representation, we discussed the Bethe ansatz of 
both theories in a unified manner. By this, in the six-vertex model at roots of unity, we obtained 
the fusion relation of -matrices, T^Q-relation from the TQ-relation, and further indicate the 
special nature of Q-operator in accord with the required s^-loop algebra symmetry of the six- vertex 
model. 
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